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ABSTRACT 



After several years of intensive technological develop- 
ment Virtual Observatory resources have reached a level 
of maturity sufficient for their routine scientific exploita- 
tion. The Virtual Observatory is starting to be used by 
astronomers in a transparent way. In this article 1 will re- 
view several research projects making use of the VO at 
different levels of importance. 1 will present two projects 
going further than data mining: (1) studies of environ- 
mental effects on galaxy evolution, where VO resources 
and services are used in connection with dedicated ob- 
servations using a large telescope and numerical simula- 
tions, and (2) a study of optical and near-infrared colours 
of nearby galaxies complemented by the spectroscopic 
data. 

Key words: Virtual Observatory; Galaxies; Stellar Popu- 
lations; Evolution of Galaxies; Photometry of Galaxies. 



1. INTRODUCTION 



What is the Virtual Observatory? The Virtual Observa- 
tory is a realisation of the e-Science concept in astron- 
omy; it is a powerful virtual environment aimed at facili- 
tating astronomical research and increasing scientific out- 
put of astronomical data. It is formed by data archives and 
software tools interoperating using a set of peer-reviewed 
standards and technologies developed and maintained by 
the International Virtual Observatory Alliance (IVOA). 

What does this really mean? Being naive, the increase 
of the scientific output of the data means that each Giga- 
byte of data coming from a given instrument will produce 
larger number of scientific results, i.e. papers or confer- 
ence presentations, exactly like uploading a research pa- 
per to the preprint server increases significantly its scien- 
tific impact. 
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Figure 1. Similarities between the concepts of WWW 
(left) and VO (right). 



2. SIMILARITIES BETWEEN VO AND WWW 



Virtual Observatory is sometimes referred as a World 
Wide Web for astronomers. Indeed, there are numerous 
remarkable similarities between the concepts of WWW 
and VO(Fig[T]i: 

• IVOA plays a similar role for the VO as W3C does 
for WWW: these are administrative organisms re- 
sponsible for the definition of interoperability stan- 
dards. As examples, we can consider specifica- 
tions of HTML/XHTML developed by W3C and 
VOTable developed by IVOA 

• Resources are the inalienable parts of both concepts. 
In case of WWW these comprise: (1) web-sites, (2) 
portals and directories, (3) web-services. In the VO 
we have: (1) data archives, (2) catalogue access ser- 
vices (e.g. SDSS CASJobs: http: //cas.sdss.org/| l, (3) 
astronomy-oriented web-services 
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• Tools represent another cornerstone of both, VO and 
WWW: (1) in WWW we deal with web-browsers 
(e.g. FiREFOx, Internet Explorer, Safari) 
while in the VO we have data browser and discov- 
ery tools, such as ASTROGRID VO DESKTOP, CDS 

ALADIN, NVO DATASCOPE; (2) advanced users of- 
ten deal with command-line tools to access the re- 
sources such as CURL or WGET for WWW, and, sim- 
ilarly, VO clients based on access libraries such as 
ASTRO-RUNTIME; (3) finally, there are specialised 
clients using WWW/VO protocols as infrastruc- 
ture and/or data transport, for instance, PICASSA or 
GOOGLE EARTH with their "analogues" in the VO 
world such as VisIVO 



3. VO SCIENCE AND TECHNOLOGY NOW 



In this section I would like to briefly mention the existing 
accomplishments of the Virtual Observatory. 

On the side of IVOA we have a comprehensive set of 
standards including: data formats (VOTable), VO re- 
source description (Resource Metadata), data model for 
ID spectra (Spectrum Data Model) and much more com- 
plex and general Characterisation Data Model, Astro- 
nomical Data Query Language, protocols to access im- 
ages and spectra, application messaging protocol allow- 
ing different VO tools to talk to each other, authorisation 
and authentication mechanisms, and others. Many more 
standards are still at different phases of development. 
Now it became possible to handle even very complex as- 
tronomical datase ts in the Virtual Observa t ory, suc h as 
3D spectroscopy dChihngarian et all l2006l l2008ah and 
results of N-body simulations. 

In the meantime, application developers have created an 
impressive set of VO-enabled tools from those of general 
interest to very specialised applications. Many of them 
were presented in the review by M. Allen (this confer- 
ence). 

Data and service providers have contributed to the Vir- 
tual Observatory by providing access to numerous data 
collections and archives at wavelength domains from 
gamma-ray to radio. First services to access theo- 
retical models (e.g. theoretical spectra of stellar at- 
mospheres in Spanish- VO or PEGASE.2 / PEGASE.HR 
synthetic stellar populations jFioc & Rocca-Volmerangel 
ll997tlLe Borgne et al.Ll2004l) in VO-France, access to the 
results of cosmological simulations in Italian VO) started 
to appear recently. We should also mention first proto- 
types of data analysis services and value-added services 
associated with data access services, such as modelling 
of the spectrophotometr ic properties of mergin g galaxies 
in the GalMer database dDi Matteo et aUl20()8h . 

The Virtual Observatory has been used for astronomical 
research for almost 5 years. The first VO-science result 
was the discov e ry o f optically faint obscured quasars by 
iPadovani et alj (12004) . This was an example of a multi- 



wavelength study carried out entirely inside the VO in- 
frastructure. Three years la ter, the VO studies o f ob- 
scured AGN were continued ([Richards et aUl2007l) . 

A number of refereed papers were published by 
the Spanish VO project members, presenting dis- 
coveries of unique objects done with the VO tools 
(ICaballero & SolanoL l2007l l2008t ICabaUero & Diiiii 
L2008|)- 

A paper byl Bayo et all (|2008|) presenting VO SED Ana- 
lyzer is of particular interest, as the first refereed paper 
presenting a "virtual instrument", i.e. a service in the VO 
aimed at data analysis, as well as its application to a par- 
ticular research project. 

Many other studies made use of the VO tools and infras- 
tructure combining them with proprietary data access an d 
analysis. For example, in lChiUngarian & MamonI ( [2008h . 
authors used the VO data discovery and access mecha- 
nisms to collect all existing data on a newly discovered 
object. This example demonstrates how difficult may be 
to define the concept of a VO study or VO-enabled study. 

Nearly all research projects mentioned above used Virtual 
Observatory to do data discovery, data access, and data 
mining. Therefore, the principal question I would like 
to address is: Is it already possible to go beyond data 
mining? The answer is: Yes, it is. 



4. VO SCIENCE BEYOND DATA MINING 

In this section I describe two VO-enabled research 
projects making heavy use of VO technologies beyond 
data mining. The VO is used in connection with ded- 
icated observations and numerical simulations demon- 
strating the proof-of-concept for such complex studies. 



4.1. Compact Elliptical and Tidally Stripped Galax- 
ies 

This study was inspired by the serendipitous discov- 
ery of a very rare compact elliptical (cE) galaxy in 
the central part of the nea rby galaxy cluster Abell 496 
dChilingarian et alll2007al) which became the 6th known 
galaxy of this class in the Universe. 

Compact elliptical galaxies have similar luminosities 
(—18 < Mb < —15 mag) and stellar masses to dwarf 
ellipticals, but 10 smaller effective radii (rg ~ 100 — 
200 pc) and, therefore, 100 times higher surface bright- 
ness and 1000 times larger stellar density per unit of vol- 
ume. The prototype of the cE type is Messier 32, a satel- 
lite of the Andromeda galaxy. All known cEs reside in the 
vicinities of larger stellar systems and/or in the innermost 
regions of galaxy clusters. 

Compact ellipticals are thought to be tidally stripped 
intermediate-luminosity galaxies dNieto & Prugniell 
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119871: iBekki et all 120011: iGrahamL l2002h . additional 
arguments for th is scena rio have been provided by 
IChihngarian et alJ (l2007ah from stellar populations. 
However, low statistics did not allow us to uniquely 
argue for this scenario of cE formation. 

Given small spatial sizes of cEs they become spatially- 
unresolved in ground-based observations for distances 
beyond ^50 Mpc. Their broadband optical colours well 
resemble K-type Galactic stars giving them little chances 
to be included into the samples of large spectroscopic sur- 
veys such as SDSS. 

The key advantage here is provided by the Hubble 
Space Telescope, which can efficiently resolve these lit- 
tle galaxies up-to a distance of 200 Mpc, allowing us to 
study their structural properties. We realised this in the 
course of our study of the Abell 496 cluster of galax- 
ies (Chilingarian et al., 2008b), and decided to search for 
cE galaxies using the power of the Virtual Observatory 
to study the role of tidal stripping in the galaxy evolu- 
tion. All details about this project will be soon provided 
in Chilingarian et al. (in prep.), here we give a brief 
overview. 
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We have constructed a workflow including the following 
steps: 

1 . querying Vizier Catalogue Service to retrieve a list 
of galaxy clusters having z < 0.055; 

2. querying NED to retrieve their precise coordinates 
and values of Galactic extinction; 

3. querying fully-reduced direct images obtained with 
HST WFPC2 and ACS from the Hubble Legacy 
Archive (HLA) using Simple Image Access Proto- 
col (SlAP); 



4. running S Extractor dBertin & Arnoutsl Il996h as a 
remote tool on these images (no image download is 
required); 

5. selecting extended objects having low ellipticity, ef- 
fective radii below 0.7 kpc and B-band mean effec- 
tive surface brightness higher than 20 mag/arcsec'^; 

6. querying NED to check if there are published red- 
shifts for the selected objects and obtaining addi- 
tional photometric data 

Having applied the workflow to the entire WFPC2 data 
collection in the HLA we ended up with the archival im- 
ages of 63 clusters with several dozens candidate cE and 
tidally stripped galaxies in 30 of them. We found a large 
number of objects in the scarcely populated region of the 
Re vs (/i)e and M vs (/x)e diagrams reflecting structural 
properties of galaxies. 

In Fig |2] we present the structural properties of newly 
discovered cE and tidally strippe d galaxy candi- 
dates in comparison with dwarf (iBinggeh & JerierJ, 



Figure 2. Structural properties of early-type galaxies. 
See sources of data in the text. New cE and tidally 
stripped galaxy candidates discovered using the VO are 
shown by red crosses on-top of green circles for objects 
with spectroscopic confirmation. 
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Gehaetal., 2003; van Zee et am2004l: iD e Riicke et al, 
20051: iFerrarese et al.L 12006 ; Chilingarian et al., 2008b), 



intermediate-l uminositv dBender et al.- ,1992). giant 
(Bender et a ll [T9^ ICaon et al.L 119931: b'Onofrio et al.L 
Il994c faberetaLL Il997h early-type galaxies, three 
nearby compact ellipticals a nd transitional cE/UCD 
object jChiUngarian & Mamonll2008l) . 

Our workflow may have confused cE galaxies with (1) 
foreground or cluster compact star-forming galaxies; (2) 
background giant early-type galaxies hosting bright ac- 
tive nuclei (AGN); (3) background post-starburst galax- 
ies. The star-forming galaxies can be discriminated au- 
tomatically by their blue colours if multi-band data are 
available or manually by clumpy morphology. The two 
remaining cases may arise when the distance to the back- 
ground object is 2-3 times the distance to the cluster, i.e. 
<400-500Mpc(z < 0.12) in our study. Then, AGNscan 
be ruled out by checking X-ray point source catalogues, 
and the probability of having a post-starburst galaxy at 
this redshift is very low (Goto, 2007). 

The next stage of the project was to obtain high-quality 
optical spectroscopic data on some of the candidates. 

We have observed three galaxy clusters, Abell 160 
(Fig O, Abefl 189, and Abell 397 hosting 8 can- 
didate galaxies (-17.0 < Mj < -19.5 mag) 
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with t he multi-slit unit of the SCORPIO spectro- 
graph dAfanasiev & Moiseevl l2005h at the Russian 6- 
m "Bol'shoi Telescop Azimutal'nyy" (BTA) in August 
2008. We have analysed the spectra by fitting them 
against high-resolution PEGASE.HR stellar population 
models usi ng the novel NBURSTS fu ll spectral fitting 
technique (IChilingarian et all l2007blH) . and obtained 
precise radial velocities, internal velocity dispersions, 
luminosity-weighted stellar ages and metallicities. All 8 
candidate objects were confirmed to be cluster members 
having stellar populations older than 8 Gyr with metal 
abundances typically between Z^jLh and the solar value 
(Zq/5 for one object) and internal velocity dispersions 
between 50 and 100 km s^^. 

We have performed numerical simulations of tidal strip- 
ping of intermediate-luminosity early-type disc galaxies 
by the potential of the galaxy cluster including the central 
cD galaxy using the GADGET-2 code (Springel, 20051). 
The simulations suggest that the progenitor galaxies may 
entirely lose their discy components due to tidal stripping, 
while keeping the bulges although with significant stellar 
mass loss as well. The remnants of the tidal stripping for 
some initial conditions and orbital configurations well re- 
semble the observed properties of cEs, although in most 
cases they are still remain quite extended. This explains 
why cE galaxies are not very common. 

We have faced a number of issues while undertaking this 
project, most of them are infrastructural. 

• There is no real VO access to NED. We had to 
develop and use "home-made" scripts to execute 
queries. 

• SIAP interface in the HLA is not publicly an- 
nounced, although it is used internally in the project. 
The service URL has been provided privately to us 
by the HLA developers, although it was possible to 
get it by reverse engineering of the JavaScript code. 

• We had to setup a customized SExtractor service and 
spend a lot of time to fine-tune it by using some un- 
documented features. 

• It took 3 semesters to convince the TAC to approve 
our telescope proposal. 

• We were the first to use the SCORPIO multi-sht 
mode with the high-resolution grism, therefore it 
was necessary to design and develop the data reduc- 
tion pipeline. 

The main result of our study is that the class of cE galax- 
ies is converted from "unique" into "common under cer- 
tain environmental conditions". We provide evidences 
for the importance of tidal stripping of stellar discs as 
a way to create peculiar early-type galaxy population in 
the cluster centres. Now we can explain the existence of 
very strange objects such as VCC 1 199 having supersolar 
metallicity for a very low luminosity of M s « — 15 mag. 



This was the first study, where the primary step was done 
in the VO, then the discovered objects were foUowed-up 
with a large telescope and successfully reproduced by nu- 
merical simulations. 



4.2. Optical and Near-Infrared Galaxy Colours 

While preparing my presentation for this meeting, I de- 
cided to make something special. The challenge was to 
get valuable results related to studies of galaxies in one 
week using the VO starting the project from scratch. My 
main collaborator in this project was I. Zolotukhin, lo- 
cated geographically in a different place, so we had to 
work remotely and interact only online. 

We decided to study optical and near-infrared colours of 
nearby galaxies and try to connect them to their stellar 
population properties. NIR magnitudes are less sensi- 
tive to the stellar population age compared to the opti- 
cal colours, therefore they can be used as better stellar 
mass tracers (although not perfect). In addition, the ef- 
fects of extinction inside the galaxies being studied are 
less important in the NIR spectral bands. Spectroscopic 
information on stellar ages and metallicities should be- 
come additional important bricks of information. The 
catalogue will be presented in Chilingarian, Zolotukhin 
& Melchior (in prep.) 

We have used the following resources: 

• SDSS DR7 dAbazaiian et aUl2008h photometric cat- 
alogues as a source of optical magnitudes 

• SDSS DR7 spectra to get stellar population proper- 
ties by the full spectral fitting 

• UKIDSS dLawrence et all |2007|) DR4 Large Area 
Survey (LAS) catalogue as a source of NIR magni- 
tudes 

The techniques we have exploited: 

• position-based cross-match (possible to do in the 
VO) 

• stellar population modelling using PE- 
GASE.2/PEGASE.HR (possible to do in the VO) 

• NBURSTS full spectral fitting technique (yet as a 
stand-alone non-VO service) 

From the zoo of VO tools we selected TOPCAT/STILTS 
to join and merge large tables, scripts-based access to 
SDSS, and ASTROGRID VO DESKTOP to access UKIDSS 
catalogues and perform the cross-match. 

We have used SDSS CASJobs to select all galaxies from 
the spectroscopic sample having redshifts 0.03 < z < 
0.3 in the SDSS stripes 9 to 16 which have been partially 
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Figure 3. HST WFPC2 mage of the Abell 160 cluster with three confirmed cE galaxies outlined (left). Spectra of these 
three objects obtained with the SCORPIO universal spectrograph { right) with their best-fitting templates as provided by 
the NBURSTS/mW spectral fitting technique. 



covered by the UKIDSS. This query has returned approx- 
imately 170 thousand objects. 

Then, we have cross-matched this list against the 
UKIDSS DR4 LAS catalogue with a search radius of 
5 arcsec. This step can be done either using WFCAM 
Archive or through the ASTROGRID VO DESKTOP appli- 
cation using the multi-cone search interface. Notice, that 
the access to UKIDSS DR4 was restricted, therefore to 
query it in an automatic way it was necessary to use the 
authorisation/authentication mechanisms provided by the 
VO. 

Then we have computed and applied iiT-corrections by 
fitting optical-NIR spectral energy distributions (SEDs) 
against PEGASE.2 stellar population models to get rest- 
frame magnitudes. 

At the final step, we have processed all selected SDSS 
DR7 spectra using the NBURSTS full spectral fitting tech- 
nique in order to estimate velocity dispersions, ages, and 
metallicities of all galaxies in 3-arcsec wide apertures. 

The non-trivial problem is the homogenisation of the 
photometric data. Firstly, we tried to use the Petrosian 
magnitudes provided in both catalogues, but we quickly 
realised that due to very different sensitivity of the two 
surveys the Petrosian radii used to measure the magni- 
tudes may be very different and the difference is corre- 
lated with the galaxy colours. Finally, we decided to 
deal with fluxes in the 3 arcsec wide apertures, which 
are provided directly by SDSS and easy to compute for 
UKIDSS using the three provided aperture magnitudes. 
These magnitudes may not reflect the real total colours of 
galaxies, since our fixed aperture corresponds to differ- 
ent spatial sizes at different redshifts, but, for explaining 
the colour properties using additional spectroscopic in- 
formation obtained by the SDSS in the same apertures, 
this approach is preferable. 



Another important problem is the iiT-correction, or de- 
pendence of colours on the redshift due to the fact that 
the filter transmission curves effectively contain different 
regions of galaxies' SEDs at different redshifts. There are 
several exist ing prescriptions for the computation of K- 
corre ctions dMannucci et al.L 120011 : iBlanton & Roweisl 
|2007|) . however, they provide controversial information 
for the NIR spectral bands. Therefore, we have decided 
to deduce X-corrections from the multi-wavelength SED 
fitting of the data against the PEGASE.2 SSP models 
redshifted according to the spectroscopic information 
and varying the effects of internal dust extinction. The 
behaviour of ii'-corrections in the optical bandpasses 
obtained in this fashion well resembles the results of 
IBlanton & Roweisl (|2007|) . At the same time, our NIR K- 
corrections turned to be very different, ho wever, well co r- 
responding to those presented in iMannucci etal.1 (1200 ll) . 

Since we had optical spectra available from SDSS for all 
the galaxies in our sample, we computed the actual val- 
ues of flux differences by integrating the spectra in the 
rest-frame and redshifted filter bandpasses unless they 
moved out of the spectral coverage. Given the wave- 
length range of SDSS spectra we were able to compute 
the "true" ii'-corrections for the g band for any redshift, 
for the r band for objects having z < 0.25, and for the 
i band till z < 0.08. These "true" values perfectly agree 
with the values derived from our SED fitting with a typ- 
ical errors not exceeding 0.05 mag. Therefore, we con- 
clude that our prescriptions for the computation of K- 
corrections have reasonable quality for the studies of op- 
tical/NIR galaxy colours. The full discussion related to 
the computation and applications of iC-corrections for 
nearby galaxies will be provided in Chilingarian, Mel- 
chior & Zolotukhin (in prep.) 

We have also faced a number of technical issues: 
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is populated by old galaxies, while in the blue cloud there 
is a significant age gradient. 

There is a 3-magnitude long high-luminosity tail of the 
red sequence clearly seen on the plots, providing an im- 
mediate evidence that such galaxies cannot be formed by 
equal-mass mergers of objects from the blue cloud. 

Another feature seen in the g — H \s Mh plot is a pop- 
ulation of young and intermediate-age objects overlap- 
ping the red sequence and sometimes being as much as 
0.7 mag redder in the g — H colour These are probably 
dusty galaxies with active star formation, so the superpo- 
sition of young and all populations, plus dust attenuation 
creates such an appearance. 

This project is still very far from being finished. We fore- 
see to add GALEX ultraviolet data and fit SDSS spectra 
together with photometric data points with more realis- 
tic galaxy models than simple SSPs, for example, to in- 
clude two star formation episodes with different dust at- 
tenuation for them. We also plan to study the distribution 
of emission line strengths since we are able to precisely 
model the stellar population, making feasible studies of 
even very faint emission lines. 



5. SUMMARY 



Figure 4. Colour-magnitude diagrams for nearby galax- 
ies: optical and NIR colours. All magnitudes are 
extinction- and K-corrected. The top and bottom pan- 
els present g — r vs Mh and g — H vs Mh respectively. 
The SSP -equivalent ages are colour coded, values from 
violet to red correspond to 1 to 15 Gyr. 

• SDSS DR7 is not accessible via VO protocols, there- 
fore we used its own CASJobs portal. 

• There were numerous problems accessing UKIDSS 
DR4 through the VO due to bugs in the implementa- 
tion of services and access interfaces. However, all 
these questions have been solved very efficiently by 
the UKIDSS and Asti-ogrid teams. 

• Need to download, upload, merge, and convert 
lengthy tables. 

• Still there is no way to perform the cross- 
match against user-uploaded table using the ADQL 
queries: these mechanisms are still to be imple- 
mented 

In Fig m we present the colour-magnitude plots for the 
galaxies from our sample. We use the rest-frame H mag- 
nitude which is much less sensitive to the effects of stellar 
age than optical bands. Spectroscopic ages obtained from 
the full spectral fitting are colour- coded. It is immedi- 
ately evident that the red sequence (IStrateva et al.L 1200 lb 



The Virtual Observatory is already at the production 
level. Scientists not associated directly to VO projects 
are trying to feel the way. The scientific results already 
obtained are impressive and very important as a proof-of- 
concept. The advantages of the VO approach are clear: 
one can transparently access and process enormous vol- 
umes of data from different sources. But, of course, the 
VO should not be considered as replacement for scientists 
- it is just a tool to help them. 

In my opinion, the major problem for a scientist in the VO 
now is very little, but numerous and, therefore, annoying 
infrastructural faults: all the individual bricks exist, but 
them together still requires a lot of efforts. 
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